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We demonstrate that hadron production viewed as formation of specific 
stochastic network can explain in natural way the power-law distributions 
of transverse mass spectra of pions found recently, which seem to substitute 
Oh, the expected Boltzmann statistical factor. 
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Recently it has been pointed out [J that properly normalized transverse 
mass spectra of vr'' mesons for niT > 1 GeV/c^ obey specific m-r scaling, 
namely they follow the universal power-law curve, 



dN 



rn^drriT \ A J 



(1) 



(which after integration over transverse momenta results in similar power 
law in masses and, according to T, describes well the yields of neutral 
mesons from r] to T, i.e., for m ~ 0.5 — 10 GeV/c^) with P ~ 8 — 10, de- 
pending on energy and type of particles. This has been regarded as purely 
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phenomenological observation of the apparent violation of the usual Boltz- 
mann behaviour, 

dN f mT\ 

- CBG • exp — , (2) 



in the region of high niT- In [2] we have noticed that this finding can be 
regarded also as confirmation that in descibing hadronization process one 
should apply not the usual Boltzmann-Gibbs (BG) statistics but rather its 
nonextensive generalization, for example in the form of the Tsallis statistics 
[H], in which eq. ((2)) is replaced by 



dN ' 



mj^dniT 



(3) 



Obviously, for large values of rriT (where dependence to scale A can be ne- 
glected) eq. becomes eq. (0). However, contrary to eq. ^ the above 
distribution offers simple interpretation of the fitted parameters, q and A. 
In particular, as demonstrated in |^ ^, g is entirely given by fluctuations 
of parameter 1/A in the usual BG approach, i.e., in eq. © ^ with nonex- 
tensivity parameter q = 1 + 1/-P. Notice that for q ^ 1 Tsallis statistics 
becomes the usual BG one and eq. © becomes eq. (HJ. 



This, however, does not solve the puzzle noticed by [J, namely what 
is dynamical origin leading to the apparent scale-free character of the ob- 
served rriT spectra seen in eq. © ^. We would like to propose here one 
possible scenario resulting in eq. (^. To this aim let us first mention that 
Tsallis distribution Q describes also very well the formation of the so called 
complex free networks 6 (if one replace uit by the number of links fe)^. 
Inspired by this observation we shall investigate in what follows the pos- 
sibility that power law seen in eq. (jSJ signals that hadronization can be 
viewed (at least from the limited perspective considered here) as process of 
formation of some specific network taking place in the environment of glu- 
ons and quark-antiquark pairs (qq), process in which their original actual 
energy-momentum distributions would be of second importance in compar- 
ison to the fact that, because of their mutual interactions, they connect to 

^ We refer interested reader to |1] 01 for details and further references. 
Notice that usually A = T, i.e., it is regarded as the "temperature" of the hadronizing 
system treated as a kind of "heat bath". Fluctuations could arise, for example, from 
the fact that, in cases considered here, such heat bath is in obvious way finite 

^ Notice that scale parameter A in eq. Q can be absorbed in the normalization 
constant, as it was done in 

* Actually, in the so called escort probability distribution was used, which results 

9 1 

in (...) 1-9 instead of (...) i-? as here. This has no consequences in what concerns 
our work as both q can be simply translated to each other. 
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each other and that this process of connection has its distinctive dynamical 
consequences. 

The proposed Une or reasoning is following. Suppose that we start with 
some initial state consisting of number no of already existed {qq) pairs, 
which we shall consider as equivalent to vertices in the usual network. We 
shall now add to them, in each consecutive time step, another vertex (i.e., a 
new {qq) pair), which can have fco possible connections to the old state. To 
be more specific, we regard our quarks as dressed by interaction with sur- 
rounding gluons and therefore somehow " excited" . Each quarks is supposed 
to interact with k other quarks, what in network terminology would mean 
that each quark has k links. We shall now assume that the "excitation" of 
quark mentioned above is proportional to k. We expect also that the num- 
ber of links k should be proportional to the number of gluons participating 
in such "excitation", i.e., existing in the vicinity of a given quark. In this 
case the natural consequence would be that the chances to interact with a 
given quark grow with the number of links k attached to it and are equal to 



i.e., that new links will be preferentially attached to quarks with already 
large values of k. This corresponds to building up of the so called preferential 
network, which evolves due to the occurence of new [qq) pairs from decaying 
gluons. Because of this the number of links is here twice the number of links 
in the usual network (cf. [Tj). After time t one has therefore 



links. This leads to the following growth equation for the i— th object (the 
i— th [qq) pair): 



where (5 = 4. Solving this equation for the initial condition stating that the 
i— th object appears in time ti with the number of links equal to k^ one gets 
that 




(4) 




(5) 




(7) 



Notice that probability of forming ki{t) < k links is then given by 
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Assuming now uniform probability distribution of occurence of objects (qq), 
the probability of adding to our system a new such object in the unit of time 
is given by 

m) = (9) 

It means therefore that, because of eq. ©, 

Pihit) < ^) = ^ (^^ > ^) = 1 - ^ (^^ < ^ • = 1 - (y )' (10) 
and probability of forming object with k links can be written as 

P(,),™<'')^_|_., („, 

i.e., the resulting distribution of the number of links is independent of time 
and has the characteristic power-like form: 

P{k) ~ k-^ (12) 

with 7 = 1 + 5 (for Q 7 = 5). For the general case of ^ 

J2ki = 2t{2ko-l) (13) 

one has 6 = = 4 — i.e., the same eq. (|T^ but with 7 = 5—^^. 

The crucial point is now our conjecture that the number of links k de- 
termines the transverse mass of emitted particle rriT- Assuming that 

rriT /c° (14) 

one gets immediately that 

_7 l-g _ 7 — 1 

P[mT) ~ mj," ■ = where /3 = 1 H . (15) 

This completes our derivation of where now P = [5. It arises because of 
the conjecture (fTH) connecting of the actual value of ttt-t with the number 
of links in some directional network growing process defined above. The 

^ It accounts for the fact that new vertex [qq] arises from the gluon — > qq process, i.e., 
one gluon disapears and this corresponds to diminishing the number of hnks by one 
for each two vertices. 

^ Notice that in the approach to networks using TsaUis statistics |5] one gets P{k) ~ 
, i.e., 7 = 5 corresponds to g = 1.25. 
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first parameter here is /cqi which denotes the number of links with which 
the new vertex (here the {qq) pair) will join the aheady existing network of 
such vertices. However, as one can see, with increasing fco one quickly ob- 
tains asymptotic value of 7 = 5. The other parameter, a, describes in what 
way the new links can be chosen in momentum spce in respect to allowed 
directions. For example, for totally random distribution one should choose, 
analogously to deterministic diffusion process, a = ^ ^. Other values of 
parameter a would correspond to other special diffusion cases (see below). 
In case of a = | one has /3 = 9 for eq. © and /3 = 9 — ^ for eq. ((T^ , which 
numerically gives us values of P equal to 7.66, 8.0 and 8.5 for, respectively, 
/cq = 3, 4 and 8. As one can see, they agree with those quoted by P (and 
listed before). 

It is worth to stress that proposed approach leads naturally to distinct 
possible behaviours of transverse mass distributions (all others will be their 
suitable combination): 

(a) The power-like distribution of the type given by eq. H12|) (which was 
our main motivation). It corresponds to the case of large excitations 
for which probability of connecting a new quark to the one already ex- 
isting in the system depends on the number of the actual connections 
realised so far. Large number of connection results then in large ex- 
citation, what means that one has large emission of gluons and what, 
finally, enhances chances of connection to such a quark. 

(6) If one assumes instead that the new quarks attaches itself to the al- 
ready existing one with equal probability (such would be situation 
for small excitations, i.e., for small px) one gets instead exponential 
distribution of links, 

P(A;)~exp(^-A^, (16) 

which results in different distributions of rriT depending on the type 
of diffusion given by parameter a in eq. p4j) and ranging from 

P(mT)~exp(^-^) for a = \ (17) 

when the full fledged diffusion is allowed to 

PirriT) ~ exp — for a = 1. (18) 

V (m-r)/ 



^ Deterministic difltusion concept is widely discused in Notice that lack of such 
diffusion in momentum space would mean that interactions are highly correlated, its 
presence indicates chaotic evolution in momentum space. 
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when there is no diffusion. This would be the case of quarks located 
on the periphery of the region of hadronization in which case they 
could interact only in interior quarks and in such case my ~ k. 

To summarize - prompted by the observation of the power-like behaviour 
of the transverse mass spectra reported by T we have considered here pos- 
sibility that they can be reflection not so much of any special kind of equi- 
librium (resulting in some specific statistics) but rather of the formation 
process, which follows (directed) free networks formation pattern discussed 
widely in the literature and found in many branches of sciences E| • Iden- 
tify vertices in such network as {qq) pairs and gluons as links we were able 
to derive eq. obtained in ^ by analysis of some experimental data, 
when we assume that observed rriT reflects somehow the number of links in 
such network (eq. (|14() ). The resulting power-like distribution is rather uni- 
versal (see, for example ,9^ where similar power-like distributions observed 
in other branches of high energy phenomenology have been atributed to the 
apparent self organizing character of the corresponding processes) and its 
power index depends mainly on the type of the allowed diffusion process as 
given by parameter a. The type of this deterministic diffusion process in the 
momentum space which is allowed in a given experimental scenario leads 
then to three kinds of distinct characteristic spectra of tut out of which all 
observed spectra could be composed. 
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